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We investigate the hydrodynamics of two identical hypersonic stellar winds 
in a binary system. The interaction of these winds manifests itself in the form of 
two shocks and a contact surface between them. We neglect the binary rotation 
and assume that the gas flow ahead of the shocks is spherically symmetrical. In 
this case the contact surface that separates the gas emanated from the different 
stars coincides with the midplane of the binary components. In the shock the 
gas is heated and flows away nearly along the contact surface. We find the shock 
shape and the hot gas parameters in the shock layer between the shock and the 
^ ■ contact surface. 

Subject headings: radiation mechanisms: thermal — plasmas — X-rays: stars 
radiative transfer — stars: neutron 

1. Introduction 

Stars of various types drive stellar winds with widely differing characteristics. In binary 
systems that consist of such stars the winds flowing out of the stars may strike each other. 
Collision of stellar winds may result in an observational appearance if the winds are rather 
powerful. Among all stars, massive OB and Wolf-Rayet (WR) stars have the strongest 
winds. The mass-loss rate M is as high as ~ 10~ 6 — M Q yr _1 for O and B stars or 
even higher (~ 10~ 5 — 10~ 4 M Q yr _1 ) for WR stars. The terminal velocity of the matter 
outflow is V°° ~ (1 — 5) x 10 3 km s _1 for OB and WR stars. The kinetic energy carried 
away by the winds is ~ 10 35 — 10 38 ergs s _1 . In the region where the two winds collide 
some part of this energy may be transformed to other forms of energy (for example, the 
thermal energy of hot gas or the energy of ultra- relativistic electrons), and then, radiated 
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at different frequencies. At present a variety of phenomena related to the collision of the 
stellar winds are observed for WR+OB and OB+OB binaries (for a review, see Williams 
1999; Cherepashchuk 2000; Corcoran et al. 2005). In particular, it is found that the X- 
ray luminosities of massive binaries are significantly higher (by a factor of ~ 10 2 for some 
systems) than the total X-ray luminosities of their single counterparts in a statistical sense 
(Pollock 1987; Chlebowski & Garmany 1991). The analysis of the energy spectra and the 
time variability of the X-ray emission indicates that the X-ray excess is emitted from the 
region where the two winds collide (Moffat et al. 1982; Chlebowski 1989; Williams et al. 
1990; Corcoran et al. 1996; Corcoran 2003). The observational data on the X-ray emission 
from massive binaries are well consistent with the model where strong shocks form because 
of the wind-wind collision (Prilutskii and Usov 1975, 1976; Cherepashchuk 1976; Usov 1990, 
1992, 1995; Stevens, Blondin, & Pollock 1992; Antokhin, Owocki, & Brown 2004). In this 
model the outflowing gas is heated behind of the shocks to temperatures of ~ 10 6 — 10 7 K 
and radiate X-ray emission that is observed as the X-ray excess. Besides, the shock waves 
in massive binaries are able to accelerate electrons to ultrarelativistic energies, and these 
electrons can generate non-thermal radio emission via the synchrotron mechanism in the 
vicinity of the shocks (Williams et al. 1990; Eichler & Usov 1993; Williams, van der Hucht, 
& Spoelstra 1994; White & Becker 1995; Dougherty et al. 1996; Jardine, Allen, & Pollock 
1996; Chapman et al. 1999; Setia Gunawan et al. 2000; Monnier et al. 2002; Pittard et 
al. 2006). For several nearby WR+OB binaries (WR140, WR146 and WR147), the radio 
sources are spatially resolved, and the radio emission model based on the wind-wind collision 
is confirmed (Moran et. al. 1989; Churchwell et al. 1992; Niemela et al. 1998; Dougherty, 
Williams, & Pollacco 2000) Dust formation may also occur in the wind collision region (Usov 
1991) and is observed for several massive binaries as an IR excess (Williams & van der Hucht 
1992; Marchenko, Moffat, & Grosdidier 1999; Monnier, Tuthill, & Danchi 2002; Williams et 
al. 2003). 

The hydrodynamics of colliding winds has been studied both analytically (Galeev, Pi- 
lyugin, & Usov 1989; Bairamov, Pilyugin, & Usov 1990; Usov 1992, 1995) and numerically 
(Luo, McCray, & Mac Low 1990; Myasnikov & Zhekov 1991; Stevens, Blondin, & Pollock 
1992; Pittard & Stevens 1997; Zhekov & Skinner 2000; Henley, Stevens, & Pittard 2003). In 
the analytical studies it was assumed that the wind of star 1 (the primary) is much stronger 
than that of star 2 (the secondary). In this case the wind-wind collision zone is near the 
secondary on the side facing the primary, and its characteristic size is much smaller than 
the distance between the components of the binary. The undisturbed gas stream from star 
1 ahead of the shock was considered as plane-parallel. The wind from star 2 was assumed to 
be spherically symmetrical. It is a reasonable approximation, for instance, to the WR+OB 
binaries where the wind momentum of the OB star is ~ 10 — 10 3 times smaller than that of 
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the WR star, i.e., rj = MobV^/MwrVwr ~ 10~ 3 — 10 _1 <C 1. However, this approximation 
is too rough for binaries with nearly the same stellar winds (77 ~ 1). In this paper we consider 
the wind-wind collision in the case when the colliding winds are identical in strength (77 = 1). 

The paper is organized as follows. In § 2 we discuss the properties of the undisturbed 
gas that flows away from the star vicinity. In § 3 we formulate the problem of collision of two 
identical hypersonic steady winds in a binary system and outline both the hydrodynamic 
equations, which describe the motion of the hot gas behind the shock, and boundary condi- 
tions on them. In § 4 we solve the set of hydrodynamic equations and boundary conditions 
and find the structure of the shock and the parameters of the hot gas in the shock layer. 
Finally, in § 5 we discuss our results and some potential astrophysical applications. 



We consider the hydrodynamics of two identical hypersonic steady winds in a binary 
system. We assume that the gas flow from the binary components is spherically symmetrical, 
nonviscous and nonheat conductive. The energy losses of the gas by radiation are neglected, 
i.e., the gas flow is adiabatic. We also neglect the rotation of the binary and its components. 

The interaction of two colliding stellar winds will manifest itself in the form of two 
shocks separated by a contact surface (see Fig. 1). The winds from the binary components 
flow radially out to the shocks. In the shock the gas is heated to the postshock temperature 
of ~ 10 7 Vs K, where V$ is the preshock wind velocity perpendicular to the shock in units of 
10 8 cm s -1 . Behind the shock the hot gas flows away from the binary system nearly along the 
contact surface (e.g., Usov 1992 and below). The contact surface separates the hot gas that 
has been emanated from the different components. For an arbitrary value of r\ the contact 
surface has a rather complex shape (e.g., Canto et al. 2006). In our case of two equal winds 
(rj = 1), there is a reflection symmetry about the midplane of the two stars, and the contact 
surface coincides with the midplane (see Fig. 1). 



We start from the set of equations, which describes the gas flow in our approximation, 



2. Formulation of the problem 



2.1. Basic equations 



div (pV) = 




p(VV)V = -Vp, 
P VV(H+ |V| 2 /2) = 0, 



(2) 
(3) 
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where V is the velocity of the gas, p is its density, p is the total pressure, 

H=^ P - (4) 

is the specific enthalpy, T is the temperature, and 7 is the ratio of heat capacities at constant 
pressure and at constant volume (e.g., Chernyi 1961). Equation (1) gives conservation of 
mass, while equations (2) and (3) conserve the momentum and total energy, respectively. 

In the shock layers between the shocks and the contact surface the hot gas may be 
considered as a rarefied totally ionized plasma. For such a plasma, we have 

p=(N e + N l )kT = ^, (5) 
m p fi 

Ni = -?-7, N e = N t Z, and (6) 
m p A 3 

where N e is the density of electrons, iVj is the density of ions, k is the Boltzmann constant, 
A is the atomic weight of ion, Z is its electrical charge, p, — A/(l + Z) is the mean molecular 
weight, and m p is the proton mass. 



2.2. Boundary conditions 

The gas parameters ahead of the shock (index 1) and behind the shock (index 2) are 
related via the Rankine-Hugoniot relations 

Pl vi n) = PsVf \ Pl + Pl [vi n) f =P2 + P2[vt h2 



2 J > 

rT/(")i2 [T/( n )l 2 
V}* = V?\ H 1 + ¥L± = H 2+ [ ^±. (7) 



Indices n and r denote the normal and tangential components of the gas velocity V. The 
condition is met on the contact surface. This condition and the Rankine-Hugoniot 

relations (7) are the full set of boundary conditions for the set of equations (l)-(5) needed 
to find the parameters of hot gas in the shock layers. 



3. Stellar winds from massive stars 



The luminosity of a massive star is very high, and the radiation pressure is responsible for 
the outflow of gas and its subsequent acceleration. The gas velocity V(r) varies from almost 
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zero at the surface of the star (r = R) to some terminal value at the distance r* ~ (3 — 5)R 
from the stellar center (e.g., Barlow 1982). 

At r > r* the gas acceleration is more or less finished, and the gravity of the star may be 
neglected. In this case, from equations (l)-(4) the spherically symmetric supersonic (OTt > 1) 
flow of gas may be described by 



v_ 

V* 



9JL 



( 7 -l)07l 2 + 2 



.(7 
(7- 



i)m 2 + 2 

1)3^ + 2" 



7 + 1 

2(7-1) 
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(7 - l)9Jt 2 + 2 

(7 - + 2 



i 
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V 



( 7 - i)m 2 + 2 
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2H _ 2 
Vf~ (7-1)^ 
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7-1 



(7 - l)97l 2 + 2 
(7 - 1« + 2 



(8) 

(9) 
(10) 

(11) 
(12) 



_(7- l)9Jt 2 + 2_ 

where 071 = V/V, is the Mach number and V s is the sound speed in the outflowing gas. The 
gas parameters marked by the sign * are taken at r = r*. For winds from massive stars, we 
have V ~ 10 3 km s _1 , V s ~ 10 km s _1 , and Wl ~ 10 2 > 1. 

For (7 — l)9Jt 2 3> 1 equations (8)-(12) are simplified and yield 

7-1 



»i, 



V = K 



V 



v°°, 
1 



p* 



2H 

V? 



«1, 

2(7-1) 



(13) 

(14) 
(15) 

( 7 -1 W i-i — ^ 

From equations (15) and (16) we can see that in a supersonic flow far from the star (r > r*) 
the effects of the thermal pressure and the thermal energy are of the order of 9JT" 2 that is 
~ 10~ 4 or less for massive stars. In our study, we neglect these effects, and the gas flow 
ahead of the shocks is described by 

P = P*(j)\ P = H = 0. (17) 



(7) 
(7) 



< 1. 



V = v c 



4. Collision of two equal winds with terminal velocities 



In Section 2 it is noted that for two colliding stellar winds the gas flow has reflection 
symmetry about the contact surface that coincides with the midplane of the binary compo- 
nents (see Fig. 1). Therefore, the problem of collision of two equal winds is identical with the 
problem of collision of one of the winds with the midplane. The latter problem is considered 
in this paper. We assume that the binary is wide (D > r*), and the terminal velocity is 
reached by the wind ahead of the shock, where D is a half of the binary separation. 



In the shock layer the position of a point may be defined by the two coordinates x and 
y that are the distances from the point to the binary axis and the contact plane, respectively 
(see Fig. 2). However, it is easier to solve the set of equations (l)-(4) using the independent 
variables ip and x, where ip is the stream function defined by the equality 



where u and v are the velocity components in the x and y directions, respectively. 

In the new independent variables ip and x the set of equations (l)-(4) can be written as 



4.1. Shock layer 



dip = puxdy — pvxdx 



(18) 



dy 1 dy v 



(19) 



dip pux ) dx u 



dv dp 
dx dip 



(20) 




(21) 




(22) 



It is convenient to introduce dimensionless variables via 



x 



r 




u = 



yoo ' 



V = 



yoo ' 
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(24) 



In these variables equations (19)-(22) take the form: 

dy 1 dy v 
dip J>ux ) dx u 

dv dp 

— = -x-L , 25 

dx di) 

_/_du _dv\ dp 
^\dx V dx) dx 1 

+ -0. (27) 

Using the boundary conditions (7) at the shock (see Fig. 2) and equation (17) for the 
gas parameters near and ahead of the shock, we can get the gas parameters near and behind 
the shock (index s): 

u s = e cos a sin(o; + if) — cos(o; + ip) sin a , (28) 
v s = — [e sin a sin(a; + ip) + cos(o; + if) cos a] , (29) 
P s = i(! -e) sin 2 (a + <p) , (30) 

H s = (1-e 2 ) sin 2 (a + v?), (31) 
where e is the ratio of the gas density ahead of the shock to that behind it, 

« = ^. < 32 > 

a is the angle between the tangent to the shock wave and the binary axis, if is the angle 
between the radius-vector r from the center of the star and the binary axis, and 



\x 2 



+ (1-I/J 2 ] 1/2 (33) 



is the dimensionless distance from the stellar center to the shock wave (see Fig. 2). For 
7 = |, we have e = |. 

The equation of the shock shape 

y = y a (x) (34) 

is connected with the angle a by 

y> • = % = tan' 1 a . (35) 
dx v ' 
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From equations (17) and (18) the stream function at the shock wave may be written as 

- f x sin (a + <p) xdx 

i> s = / =2 • (36) 

JO r s 

In the dimensionless coordinates x and ip the equation of shock shape is ip — ip s (x). 
Using equations (33) and (35), equation (36) can be rewritten as 

— f x (1 + xy' s — y s ) xdx . . 

^ S = Jo [i + (rj 2 ] i/2 [^ 2 + (i-yj 2 ] 3/2 ' ( } 

Near the contact plane the velocity component perpendicular to this plane is zero, i.e., 

v = at y = . (38) 

To find the shock shape (34) and the parameters of the hot gas in the shock layer we 
solve the set of equations (4) and (24)-(27) with the boundary conditions (28)-(31) and (38) 
by the method of Chernyi (1961) in which e is considered as a small parameter. Henceforth, 
we omit the bars over the dimensionless variables. 

We seek a solution of equations (24)-(27) in the form of a series in powers of e in the 
following form: 

y = ey , u 2 = ul + eu\, v = ev , (39) 

P = Po + epi , P = + Pi , H = H + eH 1 . (40) 

Substituting these series into equations (4) and (24)- (27) and equating the terms with the 
same power of e, we obtain the following equations for determining the terms of the series: 

w=°> (41) 

v = u dyo dyo = i ^ 

dx ' dip Pox(ul + eu\) 1 / 2 ' 
#0 = 7^-, H 1 = H [^->-\, (43) 



Hi = Hq 


(P±_Pi 


\Po Po 


poduj 


dpo 


2 dx 


dx 



(7 + 1) Po' 

dx dip ' 2 dx dx 

^(H 1 + ul)=0. (45) 
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We keep the first order term {eu\ ) in the second equation of (42) derived in zero approxima- 
tion over e. We do this because it has been shown that in this case the analytical results on 
the shock layer structure are more consistent with the results of numerical simulations (see 
Hayes & Probstein 1959; Lunev 1975 and below). 



Equations (41) and (45) are integrable by quadrature: 

u = uo(ip) , po = p (x) , 

27 p (x) 



H = Ho(ip) , po 
dpo(x) 



u\ = -2 f 
Jo 



yo 



1 n 

X Jo 



dx po(x,ip) 



(7 + 1) ff„W 



P oM)[ul^) + eu\M)} 1/2 
dy (x,ip) 



v = u (i>)- 



dx 



Pi 



-I 



^ dv {x,ifj) dip 

—~ +P*{x) , 

OX X 



H 1 = -u 2 1 (x,il>)+H*(il>), 



Pi = po(x,il>) 



Pi(x,tp) H^x,^) 



(46) 
(47) 
(48) 
(49) 
(50) 

(51) 
(52) 
(53) 



p (x) H (ip) 

where uo(ip), Po(x), H (ip), u*(if)), y*(x), p*(x), and H^ip) are arbitrary functions that can 
be determined from the boundary conditions at the shock wave and at the contact surface. 



Since y = and ip — at the contact surface, from equation (49) we have y m (x) = 0. 



4.1.1. Zero approximation 

In the series (39) and (40) we keep only the first terms marked by (zero approximation). 
To find the values u , Po, Ho, and ifj ,s in this approximation we can take that the shock 
wave coincides with the contact plane, and a is equal to tt/2. From equations (28), (30), 
and (31) we have the following boundary conditions: 

cos 2 ip TT 2 . . 

u ,s = sm tp , po iS = — — , H 0ja = cos (p. (54) 
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From equation (37), we obtain 

f x xdx 1 

^°' S = J, (1+^)3/2 =1 " (l + X 2)l/2- ( 55 ) 

Below, instead of ijj we use the new coordinate t which is the length along the shock wave, 
measured from the binary axis to the point where the stream line ip intersects the shock (see 
Fig. 2). 

At the shock wave the value of t in the zero approximation is equal to x, and substituting 
t for x into equation (55) we have the following connection between ip and t: 

^ = 1 ~ (l+t2)l/2 - ( 56 ) 

In the new independent variables t and x (0 < t < x), equations (46), (47), and (54) 
yield 

u (t) = sin (p(t) = {l+t2 y /2 , (57) 

cos 2 y?(t)| t=a; 1 ,. Q s 

Mx)= rg(x) = (iT^r (58) 

H (t)= cos 2 <p(t) = T ±-p, (59) 
r ^ 2 7 po(x) 2 7 1 + t 2 

Po (x ' t} = (7TT) luT) = (7TT) (TTW ' (60) 

where r s (x) = (1 + re 2 ) 1 / 2 [see equation (33)]. 

From equations (49) and (56) we have the coordinate y Q as a function of x and t: 

( 7 + l) (l + x 2 ) 2 /•« 

y ° ( "' t) = ^ ~/o 

_ ( 7 + l) (1 + x 2 ) 2 
47 a; 

where 



Jit) , (61) 



J(t) = arctant + . (62) 

1 ~\~ t 



From equations (50), (57), and (61), we obtain 



M^)- ^ ( 1+t 2)l/2 ^2 • ^ 
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Since t — x at the shock wave, from equations (39) and (61) we obtain the equation of 
the shock shape [see equation (34)] in the zero approximation (index z): 

y = y*(x) = ey (x,t)\ t=x = (64) 



Modified zero approximation 

In the modified zero approximation we take into account only the e correction (eu\) 
to u 2 . In this approximation we also have a = ir/2, and from equations (28) and (31) the 
boundary conditions are u l s = and s = for and Hi, respectively. Using these 
boundary conditions, from equations (48), (52), (58), and (60) we obtain 



ul(x,t) = -2j 



dp dx 
dx p 



7 (1+t 2 ) (1 + x 2 )' 



Hi(x, i) = —u\{x, t) . 



(66) 



From equations (49), (56), and (65) the equation of the shock shape in the modified 
zero approximation (index mz) can be written as 



(7-1) (1 + x 2 ) 



2\2 rx 



tdt 



x 



t 2 - 



27 a; 7 (1 + ^ 2 ) 2 

2(7-1) ln (l + t 2 )^ 1/2 



7 



(1 + x 2 ) 



(67) 



4-1.3. First approximation 

We calculate now the shock layer parameters where all e corrections in the series (39) 
and (40) are included (first approximation). 

Equations (28), (30), (31), and (35) yield the boundary conditions for pi, ui, and Hi 
near and behind the shock in the form: 

Pl,s{x) = " (1+^2)2 ' u hs = Q, H ltS = 0, (68) 
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Using the boundary conditions (68), from equations (51), (52), and (63) we obtain 

(7 + 1) (1 + 3x 4 ) f x t 2 J(t) dt 1 



(1 + t 2 ) 2 (1 + x 2 ) 2 
(7 + 1) (1 + 3x 4 ) 1 



- 87 x* -^-^-(TTW' ™ 

H^x.t) = -ul{x,t), (70) 

where 

G(x) = arctan 2 x — r- arctanx — r— r, (71) 

(1 + X 2 ) (1 + X 2 ) 2 

and Ui(x, t) is given by equation (65). 

Equations (53), (58), (59), (65), (69), and (70) yield 

o(xt)~ 4(1 + t2) ^ 1 + 3X ^ 1+X2)2 \C(X) C(t)} I ln (1 + t2) ] (72) 

From the first equation of (24) we have the equation of the shock shape in the first 
approximation (index /): 



4.2. The shock structure 

The shock shape in different approximations is given by equations (64), (67), and (73). 
From these equations it follows that the dimensionless distance from the shock wave to the 
contact plane at the binary axis (x = 0) is 



l-3e) 



1 - 2 



6 ^ 



1 + e 



(75) 



v.'(o) = (76) 

in the zero, modified zero, and first approximations, respectively. The detachment of the 
shock wave and the contact plane was calculated numerically by Lebedev & Savinov (1969) 
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and Savinov (1975), and for 1.2 < 7 < 1.8, i.e., for 0.09 < e < 0.29, its dependence on e at 
x = was fitted by 

y "(0) = (0.365 e + 0.035). (77) 

Here, we briefly discuss the shock shapes calculated in different approximations for the case 
of 7 = 5/3, which may have a special interest for colliding winds in massive binaries (see 
below) . 

For 7 = § and e = \ from equations (74)-(77) we have y z s {0) = 0.2, y™ z {0) ~ 0.106, 
y{(0) ~ 0.141, and y"(0) ~ 0.126. We can see that the e corrections in the series (39) and 
(40) make more precise the calculations of y s (0). 

The shock shapes in different approximations and the results of numerical simulations 
(Lebedev & Savinov 1969; Savinov 1975) are plotted in Figure 3 in the dimensionless coordi- 
nates x and y for < x < 1. The shock shape in the modified zero approximation is the most 
consistent with the results of numerical simulations fulfilled by Lebedev & Savinov (1969) 
and Savinov (1975). Figure 3 shows that at large distances (x > 0.9) from the axis the shock 
shape in the first approximation where all e corrections are included into consideration differs 
from the numerical result even more than the shock shape in the simple zero approximation 
where all e corrections are ignored, and this difference sharply increases with increase of x. 
Such a paradoxical behavior of the shock shape calculated in the first e approximation for 
7 — 1 ~ 1 is known and has been mentioned in many papers (for a review, see Lunev 1975). 
It is connected with the following. The e correction of the gas pressure (69) is negative and 
rather slowly decreases with increases of a:. At x ~ 1 the absolute value of this correction 
is comparable with the gas pressure in the zero approximation (58), and the used analytical 
method where all e corrections have to be small is not applicable. In fact, the paradoxical 
behavior of our solution at x ~ 1 is because for 7 = 5/3 the value of e = 1/4 is not really 
small enough. Therefore, to escape the difficulty at x ~ 1 and to improve the solution in 
the zero approximation it was suggested the modified zero approximation where only a part 
of e corrections are included (Freeman 1956, 1958). The accuracy of the last approximation 
for 7 = 5/3 is ~ 10 — 20% that is higher than the accuracy of the zero approximation (see 
Figure 3). Since in the method developed by Hayes & Probstein (1959) and Chernyi (1961) 
and used in our calculations the value of e is considered as a small parameter the accuracy 
of our results has to increases as 7 approaches unity. We hope to make a detail comparison 
between the analytical and numerical results for different values of 7 elsewhere. 
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5. Discussion 

We have considered in this paper the collision of two identical hypersonic steady winds 
in a binary system. We neglect the rotation of the binary and its components. This is 
expedient because the accuracy of our calculations is not higher than ~ 10% while in most 
observed massive binaries, for which we are going to use our results in future, typical stellar 
wind velocities exceed typical orbital velocities by a factor of ~ 10; thus, the orbital rotation 
doesn't distort substantially (more than ~ 10%) the gas flows out to distances equal to 
several orbital separations. We assume that the winds from the binary components are 
radial and spherically symmetrical. In this case the contact surface that separates the gas 
emanated from the different stars is known before calculations from a reflection symmetry 
about the midplane of the stars and coincides with the midplane (see Fig. 1). The two shock 
layers between the contact surface and the shock waves are the same because of the reflection 
symmetry, and therefore, it is necessary to consider only the properties of one of them. The 
problem is solved analytically, using the method in which the ratio of the gas density ahead 
of the shock to that behind it is considered as a small parameter, e < 1 (e.g., Bairamov et 
al. 1990; Usov 1992 and references therein). We have found the shape of the shock wave 
and the gas parameters in the shock layer. We have done this for 7 = 5/3 (e = 1/4) in the 
following three approximations. In the zero approximation we keep only the first terms in the 
e-series (39) and (40) for the shock layer parameters. In the modified zero approximation 
we keep additionally only the second term (euf) in the e-series of u 2 (39). In the first 
approximation all e corrections in the series (39) and (40) are included. We have shown that 
the shock shape calculated in the modified zero approximation is the most consistent with 
the results of numerical simulations (see Figure 3). The expected accuracy of the shock layer 
parameters calculated in this approximation is ~ 10 — 20%. These parameters may be used 
for calculations of emission from the shock layers, including the X-ray line profiles. 

The value of 7 = 5/3 has a special interest for our problem because it may relate to the 
hot gas in the shock layers in massive binaries if these binaries are not too wide. Indeed, 
for massive binaries the gas temperature in the shock layers is very high (~ 10 6 — 10 7 K 
or even higher) if the distance from the axis is not too large (x < 1). In the case of 
local thermodynamic equilibrium at such temperatures light elements (H, He, C, etc.) are 
practically totally ionized while atomic nuclei of more heavy elements (Ne, Mg, Si, S, and 
Fe) have a few bound electrons, and these strongly ionized, heavy ions are responsible for 
the X-ray emission lines observed from massive binaries (e.g., Henley et al. 2005). It is know 
that in WR winds hydrogen is essentially absent, and helium predominates. The relative 
abundances of neon and more heavier elements is < 3.5 x 10~ 3 (Dessart et al. 2000). Taking 
into account that these heavy elements are highly ionized in the shock layer, we can see that 
the deviation of the gas pressure from the pressure of totally ionized plasma is not larger 
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than ~ 1%. In this case, a totally ionized plasma is a good approximation for the state of 
the hot gas, and 7 is nearly 5/3. The same approximation may be applicable even better 
for the hot gas of the shock layers in OB binaries where the wind composition doesn't differ 
significantly from the solar abundances, and hydrogen predominates (Anders & Grevesse 
1989). Therefore, the equation of state with 7 = 5/3 is usually used for both numerical and 
analytical studies of colliding winds in massive binaries (e.g., Luo et al. 1990; Stevens et al. 
1992; Usov 1992). The question is there is or not local thermodynamic equilibrium in the 
shock layers. Recently, using an archived Chandra HETGS X-ray spectrum of the WR+O 
colliding wind binary 7 2 Velorum it was found evidence that the Mg XI emission originates 
from hotter gas closer to the O star than the Si XIII emission, which suggests that ionization 
of these elements may be non-equilibrium (Henley et al. 2005). But, this cannot significantly 
change the equation of state for the hot gas in the shock layer because the abundances of Mg 
and Si are very small. In another X-ray observations by Chandra presented by Pollock et al. 
(2005) for the wide Wolf-Rayet binary WR 140 there is evidence that the temperatures of 
ions and electrons in the shock layers are different. Such a difference of the ion and electron 
temperatures is expected for wide binaries. The point is that the postshock temperature for 
ions is much higher that the same for electrons (Zel'dovich & Raizer 2002). In the process of 
the gas outflow the electrons are heated by collisions with the ions, and their temperatures 
are equalizing if the gas density in the shock layer is high enough, i.e., if the binary is not 
too wide (Usov 1992). For two identical colliding winds the condition of the temperature 
equalization may be written as 



We can see that the restriction (78) on the binary separation 2D is rather hard if the wind 
velocity is large. For WR 140 the terminal velocities of the WR and O winds are high 
(~ 3000 km s _1 ), and therefore, the electron and ion temperatures are not equalized in the 
shock layers. If for a binary the condition (78) is satisfied the one-temperature equation of 
state (5) is correct for the hot gas, and 7 is nearly 5/3. 

It is known that the gas flow in the region of the wind-wind interaction region may be 
unstable, especially if the cooling timescale is not small in comparison with the flow timescale 
(e.g., Usov 1991; Stevens et al. 1992; Walder & Folini 2000). We are planning to take into 
account radiative cooling of the hot gas, and then, to study stability of our new solution for 
the shock layer. 

We thank the anonymous referees for many helpful suggestions that improved the paper 
considerably. The research was supported by the Israel Science Foundation of the Israel 
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Fig. 1. — Sketch of a colliding- wind binary system in which the components (hatched) are 
the sources of identical hypersonic winds. The hot gas is bounded by two shocks (solid lines). 
The gas emanated from the different stars is separated by a contact surface (dashed line). 
The directions of gas flow are shown by arrows. 
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Fig. 2. — The shock wave (thick solid line) with the equation y = y s (x) in the coordinates 
x and y where the axis x belong to the contact plane (y = 0), and the axis y coincides with 
the binary axis (x — 0). A-point is the center of the star, and N is the point where the radial 

— * 

current line AN intersects the shock wave. The vector k is the tangent to the shock wave 
at the point N, and t is the coordinate of the point N which is the length along the shock 
wave, measured from the binary axis to this point. 
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Fig. 3. — The shock shape in the zero (dashed curve), modified zero (dot-dashed curve), and 
first (dotted curve) approximations in the dimensionless coordinates x and y where the axis 
x belong to the contact plane, and the axis y coincides with the binary axis. The result of 
numerical simulations is represented by the solid curve. 



